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Photoabsorption cross sections at superhigh energies of real photons 

E. V. Bugaev a 

institute for Nuclear Research, Russian Academy of Sciences, 
60th October Anniversary Prospect 7a, 117312 Moscow, Russia 

The brief review of modern theoretical models describing the process of the photon absorption by nucleons 
at superhigh energies of real photons is given. The main aim of the work is an estimation of the theoretical 
uncertainty of the cross section prediction at photon energies around 10 19 — 10 20 eV. 



1. Introduction 

Photons of energy 10 19 — 10 20 eV are interesting 
in connection with the problem of the origin of 
high energy cosmic rays. 

The rather abundant fluxes of UHE photons 
are predicted in top-down models. UHE photons 
appear also as a result of interactions of UHECR 
with relic radiation. 

Search for UHE photons use the comparison 
of extensive air shower data with results of de- 
tailed simulations based on assumptions on pho- 
tonuclear cross section. 

Extrapolation of low energy data on a lp has 
been given by the parametrization [1] 



114,3 + 1,647 ln 2 [s/s ] (fib), 



(1) 



where s = 88.243GW 2 . The Zn 2 s-law in this 
formula has been chosen assuming the validity of 
the vector dominance model (VDM) and additive 
quark model (see below, eq. |Qj)). 

More recently, M. Block and F. Halzen [2] sug- 
gested the extrapolation of a lp using the analyt- 
icity arguments . They wrote, as a starting point 
[y is the photon energy in the lab frame): 



C a +dln- 



C 2 ln 2 —+ l 3( — r-\(2) 



Constants C, are constrained using the precise 
low energy fit at y/s < 2,01GeU (given by 
Damashek and Gilman [3J). Authors showed also 
that the fit with C 2 = (i.e., the In s - depen- 
dence in asymptotic) is not good (from a point of 
view of the x 2 -analysis) . 



2. Eikonal (minijet) models 

The total photoabsorption cross section in 
eikonalizcd minijet models is calculated by the 
basic formula (see, e.g., [1]) 



a lp (s) = 2 n/(e ) 



d 2 b x 



rlk 2 
1.2 I 1 



-x(s,b,fe_Lo)l 



(3) 



(A.- 



Here, x( s ) b, fc^ ) is the factorised eikonal, in 
which the probability that the photon can pro- 
duce the hadronic fluctuation (qq-pair) is re- 
moved, k±Q is the transverse momentum of the 
quarks of the pair. In general, the eikonal func- 
tion x is expressed through the parton densities 
inside of the hadronic fluctuation of the photon 
and inside the nucleon target, 



X (s,b,kl ) = A(b) / dp 



dxi J dx2 x 
x n l (xi,p 2 L ,k 2 ±0 )n p J (x 2 ,p 2 L )^-. 



(4) 



It is customary to separate the number densi- 
ties rii into two components: a nonperturbative 
(VMD) component (fej_o < k% - IGeV) and 
perturbative one. The latter corresponds to rel- 
atively high masses of the hadronic system pro- 
duced by the photon. Further, one can assume 
that, at least at not too high energies, the VMD 
component is dominant and perturbative compo- 
nent can be neglected. It means that partonic 
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densities (numbers of small- a; gluons) inside vec- 
tor mesons are larger than inside of the gg-pair 
with large k± . 

The well-known example of such an approach 
is the so-called Aspen model ("QCD-inspired 
eikonal model"), see [5] for a review. In this 
model the starting point is the eikonal %(s, b) for 
the case of even scattering hadronic amplitudes 
{hifpp + fpp)) which consists of three parts: 

X(s, b) = X qq (s, b) + Xq 9 (s, b) + Xgg(s, b), (5) 

corresponding to quark-quark, quark-gluon and 
gluon-gluon interactions. In particular, 

Xgg(s, b) = A(b, (J,gg)(Tgg(s). (6) 

If f g (x) ~ x~( 1+c ) (gluon structure function), one 
has u gg ~ s e . 

A factor A(b, u gg ) is the impact parameter dis- 
tribution. The parameters e, fj, gg , ... are deter- 
mined from experiment. The total cross section of 
proton-proton interaction, in the unitary (black- 
body) limit, is 

atot(s) «2 J [l-e- x »< a > v >]<Pb = 

= 2tt(— fin 2 - (7) 

Ugg SO 

(at asymptotic energies), A(b, [i)\b->oo ~ e~ Mfc . 

Now, using VMD and additive quark model, 
the eikonal of 7p-scattering can be written as [5] 

X VM °(s,b) = \a ag {s)A{b,^j2n gg ). (8) 
From here one has 

-TP = PIT • 2 /[I - e-* VMD ]<Pb ~ Jr> 2 ± (9) 
J So 

where P^ a A d ID is the probability of the photon- 
vector meson transition, P^ ID ~ 47ra em // / 2 . 

3. Regge models 

It is well known that in the Regge approach 
to deep inelastic scattering (DIS) the (effective) 
Pomeron intercept depends on Q 2 and xsj- The 
Q 2 -dependence can be connected with DIS dy- 
namics (e.g., with DGLAP evolution) as well 



as with the existence of an additional "hard" 
Pomeron. Correspondingly, there are two Regge- 
type parametrizations of the proton structure 
functions, which are frequently used. 

In Donnachie-Landshoff model of two 
Pomerons [6] one has, in the Regge limit, 

F 2 (x Bj ,Q 2 ) = ^Ai ( -^—) 1+€i x^ , (10) 

i 

and, in the photoproduction limit, a lp is finite. 
The example of the fit is [6| [2v = s[GeV 2 }): 

a lp = 0.283(2^) - 418 + 65.4(2z/) a0808 [ub]. (11) 

Authors warn against of using the eq. (fTU)) at ex- 
tremely small xsj where the fixed-power behavior 
will be moderated by shadowing suppression. 

In the CKMT model [7] the main statement is 
that the rescattering (absorption) corrections in 
applications of the Regge theory to DIS at HERA 
energies are not small. The Pomeron intercept 
1 + eo = 1.0808 is not the true Pomeron inter- 
cept (i.e., it does not correspond to the "bare 
Pomeron"), but rather is the effective one. The 
relative contribution of the most important ab- 
sorptive corrections depends on Q 2 . As a result, 
at large Q 2 we see the bare Pomeron (with inter- 
cept - 0.25) and at Q 2 = wc see the Pomeron 
with effective intercept 1.0808 (i.e., in the photo- 
production limit, there is no term fastly growing 
with the photon energy). 

One should mention also the Regge-type model 
of [H] , where the concept of the effective Pomeron 
intercept is used, the value of which is weakly 
dependent, in the photoproduction limit, on the 
photon energy. 

4. Colour dipole models 

Real photons and virtual photons with small 
XBj have hadronic properties. The length of 
hadronic fluctuations is large (and larger than the 
target size) at large energies, I ~ — = — for vir- 

° < G ° 5 nipXBj 

tual photons and I ~ for real photons (M qg 

is the invariant mass of the hadronic fluctuation) . 
So, the 7,'V-scattering is the two-step process. 

The basic equations of the colour dipole model 
are very simple due to the fact that colour dipoles 
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Figure 1. The photoabsorption cross section as 
a function of photon energy, for different models: 
CKMT 0, DL 0, GBW [13], Aspen [5], BH [2], 
BB p, PT 0. 



are eigenstates of interaction in QCD and the well 
known method of eigenstates [9] can be used. The 
total 7p- cross section is expressed through the 
total gluon-nucleon cross section, 



d 2 r_L* 7 (r_L)6-(r ± ). 



(12) 



Here, ^ 1 (r±) (more exactly, *t,l(' , _l, z, Q 2 )) is 
the photon-dipole wave function which depends 
on the photon virtuality, on the longitudinal mo- 
mentum fraction z carried by the quark, and on 
the dipole size r±. These variables are "frozen" 
during interaction. 

The dipole cross section, &(r±), is assumed to 
be flavour-independent and depending, except of 
rxj on s = W 2 or x' , where x' is the momen- 
tum fraction of the proton carried by the gluon 
attached to the qq-\oop. 

Perturbative QCD leads to the formula [TU] 

a(r x ,x') = yria s (Q"V 5 (x',Q 2 ), (13) 

where Q 2 is the energy scale depending on the 
dipole size, g(x',Q 2 ) is the gluon density. 

In MFCS model [IT] it is assumed, phenomeno- 
logically, that Q 2 = 4- , A = const (4 10). 



Using unitarity, as a guide, the dipole cross sec- 
tion can be written as 

a(r x , x') = 2 J d 2 b[l - e- x( - b ' E ^ r ^] = 

= 2 j d 2 6r(£; 7 ,rx,fe). (14) 

The profile factor T(E 1 ,r±,b) is smaller than 1. 
If a(r±,x') is known one can calculate T assum- 
ing some law of i-dependence of the scattering 
amplitude. 

T.Rogers and M.Strikman [12] calculated, us- 
ing MFCS model, the jp - cross section up to 
super-high energies. They unitarized cross sec- 
tion "by force", calculating profile function T, 
and, if the profile function exceeded unity, putting 
r = 1. By such a way they determined the max- 
imum rise of a lp with energy, <j 1p ~ In 3 _E 7 , for 
E 1 > !0 3 GeV. 

The dipole cross section <j(r± 7 x') in MFGS 
model rises with energy infinitely. If we suppose 
that this cross section is bounded by an energy 
independent value, as in simple saturation mod- 
els [13], the rise of a lp with energy still will take 
place, due to the photon wave function. In the 
model of [13] one assumes that 



(15) 



<T(x',r ± )=*o[l-e 4 *o<*'>], 

R 2 (x') = (x'/xo^GeV- 2 , x' - x Bj - If the 
Bjorken variable, in the photoproduction limit, 
is modified to be 



X = XB 3 (l 



4m? 
Q 2 ' 



4m| 
W 2 



(0 2 ^o), 



(16) 



it is easy to show that a lp ~ In E 1 . 

The modification of the model [13], with tak- 
ing into account the QCD evolution of the gluon 
distribution [T4] . 



a(x',r±) = ct [1 - e 



(17) 



(here the scale p 2 is assumed to have the form 
p 2 — + Pq,), gives, in asymptotics, a lp ~ 

In 1 ' 2 

Physically, the saturation model, in its phe- 
nomcnological variant, corresponds to the proton 



4 



E.V. Bugaev 



being a disc in the transverse plane with a sharp 
bordering. Saturation leads to an uniform black- 
ening of the disk with decreasing XBj (for a parent 
dipole with a fixed size r±) without changing the 
disc size. 

For a GBW model [1 3j the parameters are: 



A = 0.3 ; io = 3x 10" 



(Tq = 23m6. 



(18) 



5. Unitary bounds on s-dependence of o lp 

Using three assumptions: i) 7_/V-interaction is 
a two-step process, ii) generalized VDM (a dis- 
persion relation with variable M 2 , where M is 
the mass of the hadronic fluctuation), in) the 
hadronic interaction of the fluctuation (gg-pair) 
is a black disc interaction, one can show that the 
<7 7 * n is given by [15] 

R(M 2 )M 2 dM 2 



3tt 



{Q 2 + M 2 f 



-o- M 2jv(s), (19) 



where R 



cr(e 4 



^hadrt 



If we suppose that 



cr(e+ e — >/-t + M ) 

the hadronic cross section ctjv/ 2 jv is given by the 
black-body limit, <Jm 2 n ~ In 2 — , one obtains, in 
the photoproduction limit, 



2 s M, 
c 7 at ~ m — m 



Q 2 



in 



Q 2 



In*—. (20) 



It is the so called Gribov bound [15]. However, 
the hypothesis of the black disk interaction can- 
not be correct for gg-pairs with large mass M q q. 
Such a pair has a small rj_ ~ 1/M q q and, being 
colour neutral, interacts with the target weakly, 
a ^ r\ (it is predicted by pQCD). Due to this, 
there is the following constraint on the value of 
M [IT]: 



M 2 



^2 e ^/aln(l/x Bj ) 



(21) 



(A is the QCD scale), so, in far asymptotics one 
obtains the corrected unitary bound (assuming 



that x 



\/s, as in eq. (16)) [T7] : 



a lN ^ln 2 (^)ln^ 2 (4-)^ln 2 ^. (22) 
\s J Vs™"/ s 



trophically. At s — 10 n GeV 2 the predictions are 
in interval (0.5 -j- I.l)m6. 

It is rather difficult to predict reliably the 
asymptotic s-dependence of the a lp . Most prob- 
ably, the upper limit on the rise of the photoab- 
sorption cross section is given by the law ~ In 
as follows from the estimates of authors of [12] 
based on the colour dipole model. The corrected 
Gribov unitary bound based on generalized VDM 
and pQCD constarint gives a lv ~ ln 2 - 5 E 7 . Pre- 
dictions based on the vector meson dominance 
and additive quark model give the law ~ ln 2 E 1 . 
At last, Regge eikonal model [7] predicts, asymp- 
totically, <7 7P ~ E®' 1 . 
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6. Conclusions 



The straggling of theoretical predictions for a JP 
inl9 in20„T/ ;„ i„^„„ 
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at E~, ~ 10 19 — 10 20 eU is large, but not catas 



